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Abstract 

In this paper we describe new classes of periodic solutions 
for point vortices on a plane and a sphere. They correspond to 
similar solutions (so-called choreographies) in celestial mechan- 
ics. 



1 Equations of motion and first integrals 
for vortices on a plane 

For n point vortices with Cartesian coordinates (xi,yi) and intensi- 
ties Tj, the Hamiltonian equations of the motion are 

ViXi = ^ T i y i = -^L, l^i^n, (1) 

% OXi 

where the Hamiltonian is 

n 

ft = -j- z2^j In \n - rj\ 2 , ri = (xi,yi). (2) 

i<j 



For these equations, the Poisson bracket is {xj, yj} = F i 1 5, 
The system has three first integrals 



Q = J2^ P=J2 T ^ / = Z) r <(^ + 2/i), ( 3 ) 



i=l i=l i=l 
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which are not involutive: 

JV 

{Q,P} = ]Tr,, {P,I} = -2Q, {Q,I} = 2P. (4) 

i=i 

However, two involutive integrals always exist, for example, Q 2 + P 2 
and /. With the help of these integrals, we can reduce the system by 
two degrees of freedom. 

Thus, in the case of three vortices the reduced system has one degree 
of freedom and it is integrable (Grobli, 1877 ^3]; Poincare, 1893 |2*T]). 

The four- vortex problem is reduced to a system with two degrees of 
freedom. This system is not integrable (S. Ziglin, 1979 [2*7]). 



2 Equations of motion and first integrals 
for vortices on a sphere S 2 

For n point vortices moving on a sphere S 2 , with spherical coordi- 
nates (8i, (fj) and intensities Tj, the Hamiltonian equations can be writ- 
ten as (Bogomolov, 1977 [I]) 

9 i = {9 il H}, & = {^i,H}, z = l,...,n, (5) 

with the Poisson bracket {(fi, cos 9k} = * fc and the Hamiltonian 

R r,- 



i<k \ ' 



(6) 



Here, R is the sphere's radius, and 7^ is the angle between the vectors 
from the center of the sphere to the point vortices i and k so that 

cos 7^ = cos 9i cos 9f. + sin 9i sin ^ cos(y?j — fk)- 

Beside the Hamiltonian, the equations (0) have three independent 
non-involutive integrals 

n n n 

F\ — R. r i sin 9i cos ipi , F 2 — R Tj sin 6*j sin ipi , F3 = R T j cos 9i . 

i=l i=l i=l 

(7) 
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The vector F = ^2 ^i r i i r i are the radius- vectors of the vortices) 
with components (Fi, F 2 , F 3 ) is called the moment or center of vortic- 
ity. The first integrals commute in the following manner: {Fi,Fj} = 

= —EijkFk. Similar to the case of a plane, we can reduce the system by 

two degrees of freedom, using involutive integrals, for example, F 3 , \F\ . 

So, in the case of three vortices, we have a completely integrable 
system (P.Newton and R. Kidambi, 1998 [Ej, 2000 j2D]; Borisov and 
Lebedev, 1998 jH]). The four- vortex problem is reduced to a system 
with two degrees of freedom and, generally, is not integrable (D. A. and 
A. A. Bagrets [3]). 



3 Reduction of a system of vortices on a 
plane 

To find new periodic solutions and apply the numerical analysis, we 
will do the most symmetric reduction of the system by two degrees 
of freedom. We present here one of the possible ways to do such a 
reduction, based on new equations of motion for the mutual variables. 

Effective reduction of a four- vortex system was done by Khanin ^7] 
(for vortices of the same sign) and by Aref and Pomphrey [2] (for equal 
vortices). The generalized Jacobi reduction by one degree of freedom 
in the point vortex dynamics was described in |19j . 

In our study, we used the universal algebraic method of effective 
reduction (for an arbitrary number of vortices). It is based on mutual 
variables representation of the motion equations of a system of point 
vortices 0, PD- 
As mutual variables, we introduce the following quantities. Squared 
distances between the pairs of vortices M^- and oriented areas of trian- 
gles Ayfc, 

M i:j = (xi - Xjf + (iji - yj) 2 , A ijk = (rj - r k ) A (r k - r t ). 

These variables are due to E. Laura. 

We put Ti = Tj • = 1, P = Q = 0. Then the momentum integral I (J3J) 
can be written as 

n 

J = ^E M ^ ( 8 ) 

i<j 

where n is the number of the vortices. Writing the vortices' coordinates 
in the complex form z k = x k + iyk, we obtain the following representa- 
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tion for them: 

n 

** = £Ev^*e*«, ( 9 ) 
j=k 

where O^j is the angle between the vector from the j'th vortex to the kth 
vortex and the positive direction of Ox. 

By direct calculation we can derive the following propositions, which 
describe the behavior of the reduced system of three and four vor- 
tices 

Proposition 1 For three identical vortices, the evolution of the mutual 
distances (with fixed I = const) is described by a Hamiltonian system 
with one degree of freedom. In terms of the canonical variables (g,G), 
this system can be written as 

9=jjg, ^=-fr> H = -±\nM 12 M 13 M 23 , (10) 



where M 12 = 4(|-<?), M 13 = 8G-I + 2y/T2~J (j^-G)G cos g, 



M< 



2:i 



4(|-G) - 2VT2MI - G) Gcosg. 



Proposition 2 For four identical vortices, the evolution of the mu- 
tual distances is described by a Hamiltonian system with two degrees of 
freedom. In terms of the canonical variables (g, G, h, H), this system is 

g -dG> h ~dH> ~~dh' (U) 

U = --^\nM 12 M 13 M u M 23 M 2i M u , 

where 



M12 = I-G+2y/(I-H)(I-G) cosh, M 34 = I-G-2y/(I-H)(I-G) cosh, 
M13 = I+G+2^(I-H)Gcos(h+g), M 24 = I+G-H-2 ^(I-H)G cos(h+g), 
M u = H+2 y/(H-G)G cos g, M 23 = H-2y/(I-G)G cos g. 

Remark 1 The above canonical variables naturally follow from the 
Lie-algebraic representation of the equations of motion. 
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Figure 1. The phase portrait of the reduced system in the three- vortex 
problem, and the relative choreographies, corresponding to different phase 
orbits in the portrait. 



4 Absolute motion — quadratures and ge- 
ometric interpretation. 

According to when Mij(t) are known, one needs to find the an- 
gles 9{j(t) in order to determine the coordinates of the vortices. Clearly, 
only one of these angles is independent (let it be 6*12), the remaining 
are obtained using the relations 

da + @ik = arccos ( — — — = — j , % ^ j, k ^ i. (12) 

The evolution of is found by integrating the first order equation: 

^ = ^7 E r *+ E Tk + 1^) - ^ E r * + ^) • 

(13) 

There is an interesting geometric interpretation of the absolute mo- 
tion of the periodic solutions of the reduced systems (fTOj) . fTTf . 
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Proposition 3 Let j(t) be a periodic solution (of period T) of the re- 
duced system. Then 

1° there exists a frame of reference, rotating with the constant an- 
gular velocity Q a about the center of vorticity, where each vortex moves 
along some closed curve 

2° the rotating velocity tt a is given by (accurate to — q, p, q E 7*): 

T 

Q - = fj 12 (t)dt; (14) 
o 

3° if f2 a and £l = — are commensurable (i. e. ^ = ^, p, q G Z), 

then in the fixed frame of reference the vortices also move along closed 
curves; 

4° if any of the curves transform into each other by rotation 
about the center of vorticity by an angle, commensurable to 2n, then 
there is a (rotating) frame of reference where the corresponding vortices 
orbit one and the same curve. 

PROOF 

In the right hand side of (|13|). there is a periodic function of pe- 
riod T. We have a convergent Fourier expansion: 

2ivn 

W 12 = ^fl( n )e l ~'. (15) 

Integrating (JT5j) . and using (fT2")h we find that all the angles % depend 
on t in the following way: 

o.jii) a ; / • //,,(/). (i6) 

where tt^ = a<® + Pv e Z ' while 9ij{t) = 9ij{t + T ) are T- 

periodic functions of time. 

Substituting the above dependency into Q, we find that the posi- 
tions of the vortices on the plane are given as follows: 

Zk {t) = ^±^ + u k {t)e^\ u k (t) = u k (t + T 1 )eC, n = a<°>. (17) 

It follows that in the frame of reference, rotating about the center of 
vorticity with the angular velocity Q, all the vortices describe analytic 
closed curves, given by functions u k {t) G C. 

The proofs of 2°, 3°, and 4°, using the relations (fTKj) - (fT7]) . are ob- 
vious. ■ 
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Remark 2 Without any modification, Proposition [3] is generalized to 
the case of n vortices if we suppose that j(t) is a periodic solution of a 
reduced system with 2n — 2 degrees of freedom to the n-vortex problem. 



5 Particular solutions and stationary con- 
figurations for n equal vortices 

This section presents most well-known particular solutions for point 
vortices on a plane and on a sphere. 

The case of M 2 : 

1. collinear configurations; 

2. polygonal configurations (in particular, polygonal configurations, 
embedded into each other). An n-polygon (Thomson configura- 
tion j2H]) is stable if n ^ 7, and unstable if n > 7 (the final result 
on the stability of the heptagon was obtained recently in [TH]); 

3. non-symmetrical stationary configuration (for n ^ 8) [T3] . 

The case of § 2 : 

1. Collinear configuration. 

The roots 6i that define a collinear configuration, rotating with 
velocity Q, can be found as equilibriums of a system of N particles 
on a circle with the Hamiltonian 



This system was obtained and studied by A. V. Borisov and V. V. Kozlov 
(1999 [Bj) — a special case with SI = is called the Dyson sys- 
tem. It is proved that the system (fT%|) is not integrable, but a 
quasi-integral exists, which gives a very good approximation of 
the behavior of systems at low energy. 

2. Analogs of Thomson configurations 0. 

In these configurations, the vortices are positioned at the same 
latitude 6q in the vertices of a regular n-polygon and rotate about 
its center with the angular velocity 




k=l k=l k,i=l 




(18) 
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r(iv-i) ctg^ 

AttR 2 sin 6*0 

6 Choreographies 

Recently, a new class of periodic solutions to the classical n-body prob- 
lem has been obtained in celestial mechanics (C. Moore, R. Montgomery, 
A. Chenciner, C. Simo, and others). For example, they found a remark- 
able stable periodic orbit in the 3-body problem — namely, a figure- 
of-eight \12\ . In this solution, three bodies follow each other along the 
same curve. 

A long and ingenious analytical proof was done using the variational 
method based on the least action principle. These solutions where all 
particles move along the same curve are called choreographies jTTj. 

Choreographies can be: 

a) absolute (in the fixed frame of reference) or relative (in rotating 
frames of reference) ; 

b) simple (a single closed curve) or complex (more than one closed 
curve) 

Similar solutions exist for point vortex dynamics on a plane and on 
a sphere ^Oj- However, the least action principle cannot be applied to 
an analytical proof in this case. 

Using Proposition 01 and the known analytical solutions [15] . one 
can prove the following theorem for three-vortex systems: 

Theorem 1 If, for the motion of three vortices of equal intensity, the 
constants of motion, I and 7i, satisfy the inequalities 

-ln3 < ^-H + \nI < ln2, (19) 

then such motion is a simple relative choreography (see Fig. 11). 

Remark 3 Note also that various properties of the motion of three 
vortices are discussed in liffi and For example, the paper 

focuses on the study of stability of collinear configurations. However, 
properties of the absolute motion have not actually been studied. 

The four-vortex system has a remarkable symmetric analytical solu- 
tion — Goryachev's solution ^3] (see also PP) where the vortices form 
a parallelogram as they move on a plane. 
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Theorem 2 If, for the motion of four vortices, the vortices (of equal 
intensity) form a centro symmetrical configuration (a parallelogram), 
while the constants 7i and I satisfy the inequalities 

-In 2 < ^H + lnl < -m±#, 
6 5 

then the motion is a simple relative choreography (see Fig. HJ). 
For the complete proofs, see [TP] . 

Let us show that for the system of four identical vortices, there exist 
choreographies, which do not coincide with Goryachev's solution. These 
are non- analytical choreographies. These solutions can be obtained by 
variation of the Hamiltonian near Thomson's solution. 

Indeed, near the Thomson configuration, the Hamiltonian can be 
written in as 

n = n 2 + n*, de g ?0 3, 

where 7i 2 is a normal form 

H 2 = ±(3(x 2 + P 2 x )+2V2(y 2 + Pl 

and Ti* is the perturbation function (a term of higher order). 

On the level Ti, 2 = h 2 = const, there are precisely two nondegenerate 
periodic motions of the Hamiltonian system with the Hamiltonian 7i 2 : 

x = P x = 0, y 2 + pi = JI-h 2 , (20) 

V2 

y = P y = 0, x 2 + Pi = ^-h 2 . (21) 

According to the Lyapunov theorem, these periodic solutions are 
conserved under the perturbation 7i* and near Thomson's solution 
there are two periodic solutions of the whole system 7i 2 + 7i*: 

solution (J2T1|) gives Goryachev's solution, and 



solution (j21)l gives a non-analytical relative choreography such as 
shown in Fig. El 

Remark 4 Here a solution is called non- analytical if it cannot be ex- 
pressed explicitly in terms of quadratures. As a rule, this means that 
there are no explicit symmetries for these solutions, since all known 
symmetrical solutions can usually be found in terms of quadratures. 
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OOO 



E= E T 



E= 1.35, n = 5.945 E= 1.364 « £ e = 6.635 



Figure 2. Relative choreographies for Goryachev's solution 




1.215, Q = 5.118 



i.28. n = 



Figure 3. Relative choreographies for the solution (1211) . 

7 Absolute choreographies 

Let us discuss in more detail, whether absolute choreographies exist in 
the three- and four-vortex problems. According to the above, any rel- 
ative choreography, corresponding to a periodic solution (of period T) 
of the reduced system (see (fTU|) . (fTTJ)) closes in time mT, m G N. Dur- 
ing this interval the vortices pass through the same configurations m 
times. We denote the corresponding angular velocities of the frames of 
reference by Qm ■ 

For three- vortex choreographies (Theorem [I}, there exists a (rotat- 
ing) frame of reference, where a choreography closes in the smallest 
possible time T. For the case of four-vortex choreographies (Theo- 
rem^, this time equals 2T. We denote and the corresponding angular 
velocities by and . The angular velocities of other relative 
connected choreographies are now given by the following relations. 

For three vortices: 

n^\E) = tt?\E) + f|fio(£), where m e N, k G Z; 



here 3k and m are coprime numbers. 
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E = 0.30257 
b) 



Figure 4. The angular velocities of the frame of reference for relative chore- 
ographies of three vortices; two absolute choreographies corresponding to 
the points a and b. 



For Goryachev's solution: 



where m is odd, fcGZ, the numbers k and m are coprime, while Qo(E) = 
= j^)' Here, t ne choreography that closes in mT corresponds to the 

velocity £1$ (and the choreography that closes in 2mT corresponds to 
the velocity f4m)- 

Figs. 4 and 5 show the angular velocities in terms of the energy pa- 
rameter for three and four vortices, respectively. The roots of flm(E) = 
= correspond to the absolute choreographies shown in Figs. 4 a, b 
and 5 a, b. It can be shown that there is a countable set of energy 
values such that choreographies described in Theorems ^ and |2] are 
absolute and simple. 
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Figure 5. The angular velocities of the frame of reference for relative chore- 
ographies of four vortices; two absolute choreographies corresponding to the 
points a and b. 

8 Choreographies on § 2 

It can be shown, using the above-described methods and the reduction 
described in j7j, that there are similar relative and absolute choreogra- 
phies in the systems of three and four vortices on a sphere. Several 
examples of spherical choreographies of four vortices are given in Fig. 6 
(they are quite similar to those shown in Figs. 2 and El). 

In [22] and periodic (in a fixed frame of reference) solutions 
for different numbers of vortices on a sphere are specified, which admit 
different discrete symmetry groups. These motions are disconnected 
choreographies. Using the group-theoretic methods, one can construct 
a large number of disconnected choreographies in vortex dynamics on 
a sphere, generalizing those specified in [22] and [2H]. However, the 
choreographies that we found cannot be obtained using these methods. 
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